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We investigate quintom evolution in power-law potentials. We extract the early-time, tracker 
solutions, under the assumption of matter domination. Additionally, we derive analytical solutions 
at intermediate times, that is at low redshifts, which is the period during the transition from matter 
to dark-energy domination. A comparison with exact evolution reveals that the tracker solutions 
are valid within 98% for z > 1.5, while the intermediate ones are accurate within 98% up to z ~ 0.5. 
Using these expressions we extract two new u)-parametrizations, one in terms of the redshift and one 
in terms of the dark-energy density-parameter, and we present various quintom evolution sub-classes, 
including quintessence-like or phantom-like cases, realization of the — 1-crossing, and non-monotonic 
' lo-evolution. 
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INTRODUCTION 



O According to cosmological observations the universe is experiencing an accelerated expansion, and the transition 
to the accelerated phase has been realized in the recent cosmological past In order to explain this remarkable 
u: behavior, one can modify the theory of gravity or construct "field" models of dark energy. The field models that 
(-h | have been discussed widely in the literature consider a canonical scalar field (quintessence) [H, 0, , a phantom field, 
O that is a scalar field with a negative sign of the kinetic term 0, or the combination of quintessence and phantom in a 
unified model named quintom Q . Such a combined consideration intends to describe the crossing of the dark-energy 
equation-of-state parameter w through the phantom divide —1 which is not possible in simple quintessence or 
simple phantom scenarios where w lies always on the same side of that bound (w > — 1 for quintessence and w < — 1 
C^j \ for phantom). 

Although one can find potential- independent solutions in field dark energy models |9| , in general the cosmological 
evolution depends significantly on the potential choice. The power-law potential is a widely studied case since, amongst 



others, it allows for a theoretical justification through supersymmetric considerations [1(|. For quintessence [3. [ill, 
or phantom case [l3| it has been shown that it exhibits "tracker" behavior jthat is for a large class of initial conditions 
the cosmological evolution converges to a common solution at late times which can be extracted analytically. At 
intermediate times the dark-energy scalar-field is still small, but non-negligible, and one has to rely on perturbative 
analytical expressions fl2l . [T3 |. However, such an extension to low redshift is very useful, since observations like 
supernovae la, WMAP and SDSS ones, are related to this period [lj. Finally, an additional advantage of these models 
is that the field energy density remains small at early and intermediate times and thus the known cosmological epochs 
ON , are not disturbed. 

In this work we desire to study the quintom scenario in power-law potentials, both at high and low redshifts, and 
provide the tracker solutions and the perturbative analytical expressions respectively. Due to the complex nature of 
the model, we expect that the results will be qualitatively different from simple quintessence and simple phantom 
models, even at the tracker level. The plan of the work is as follows: In section [II] we construct the quintom paradigm 



in power-law potentials and we extract the matter-dominated solutions. In section IIIII we derive the cosmological 
solutions at intermediate times, that is when dark energy is non-negligible, but still sub-dominant comparing to 
the matter content of the universe. In section IIVI we compare our analytical expressions with the exact numerical 
evolution of the quintom scenario, and in section [V] we provide some applications of the model, extracting two new 
w-parametrizations and discussing the cosmological implications. Finally, section [VTl is devoted to the summary of 
our results. 
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II. EARLY TIME COSMOLOGICAL EVOLUTION: TRACKER SOLUTIONS 

We start by constructing the simple quintom cosmological scenario in a flat spacetime. The action of a universe 
constituted of a canonical <j> and a phantom er field, is [3]: 



S = I d^x^J—g 



(1) 



where we have set 8nG — 1. V^cp) and V a {a) are respectively the quintessence and phantom field potentials, while 
the term C m accounts for the (dark) matter content of the universe, considered as dust. Finally, although we could 
straightforwardly include baryonic matter and radiation in the model, for simplicity reasons we neglect them since we 
are interested in z < 20 era. The Friedmann equations read [7[: 

H 2 = l(pm+ P4>+ Pa), (2) 
H = -\{pm + P4> + P<f> + Pa + Pa), (3) 

and the evolution equations for the canonical and the phantom fields are: 

p^ + 3H(p^+p^ =0, (4) 
Pa + 3H(p a +p a ) =0, (5) 

where H = a/a is the Hubble parameter and a = a(t) the scale factor. In these expressions, p m is the dark matter 
density, and and p^ are respectively the pressure and density of the canonical field, while p a and p a are the 
corresponding quantities for the phantom field. They are given by: 

= ^ 2 + W) (6) 
P4> - ^ 2 -W), (7) 

Pa = -\d 2 + V a {a) (8) 
Pa = -~a 2 -K((7). (9) 
Using these expressions, we can equivalently write the evolution equation in field terms as: 

^ + 3^+^0=0 (10) 

a + 3H*-^± = 0. (11) 
oa 

The equations close by considering the evolution of the matter density, which in the case of dust reads simply 
p m + 3Hp m = 0, leading to p m = p m o/a , with p m Q its value at present. 

In quintom scenario, dark energy is attributed to the combination of <j) and a. Thus, we can define the dark energy 
density and pressure as pde = P<t> + Pa and pue = P4, + Pa- Therefore, the dark energy equation-of-state parameter 
reads 

Pde Pcj, + Pa , 10 s 
w = = — . (12) 

Pde P0 + Pa 

In the present work we study quintom evolution in power-law potentials. That is we consider: 



V^(<f>) = ^' a * (13) 
V a {p) = K a cj- a % (14) 
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where and n a are constants with units 7n 4+Q * and m 4+Q<T respectively. Note that we do not restrict the values 
of and a„ a priori. However, as we are going to see, energy positivity will force <x$ to be positive and a a to be 
negative and bounded, that is potential (fl3|) will be an inverse power-law one while ([Ft)) will be of a normal power-law 
form. 

We can express the aforementioned cosmological system using the scale factor a as the independent variable, since 
it is straightforwardly related to the redshift z which is used in observations. Generalizing [12L [l3l] we define 

P<t> +P<t> \<P _ 1 / d 4>\ 2 hK \ 
* 2(p m + p<>) SH 2 6\ a da) ' UOj 

- - Pa+Pa = -l(a^)\ (16) 



and thus (O,© can be written as: 



Therefore, we can simply write: 



and thus for the dark energy equation-of-state parameter we obtain: 

x^pm — V^((t))(l — 2x^) x a p m -V a (a)(l-2x a ) 



2{p m + p a ) 3H 2 6 V da J 



U 2 = T^r[pm + vm, (17) 

-t^Y^iPm + V^a)}. (18) 



p a = ^PHL + W (20) 
1 - x a 

^ = T—^ (21) 

x a p m -V a (a)(l~2x a ) 

Pa = : , (22) 



W ~ x^pm+V^) x„p m + V a ~) • ^ 6 > 



1 X 1 X (j 

Similarly, the first Friedman equation can be written as 



r2 _ Prn + W) + V a {<j) 



■AH 2 = ' m ' ' ay (24) 

1 X(p X (j 



Finally, using expressions p^|) - ([^)) and the field evolution equations (TrTj|) and (fTT|) become 

00 + - (5 - 3 ^- 3 ^ ) + p m + W)+K(a)l 2 + —\=° (25) 



a a" + -— (5 — 3x0 — 3xc) + 



3(1 - - jg) 

Pm + W)+K(<x) 



/[y (0)+K(a)] dU CT (a) 



= 0, 



(26) 



with prime denoting the derivative with respect to a. These equations are exact and account for the complete dynamics 
of the quintom scenario. 

In this section we desire to present the tracker solutions, that is the initial-condition-independent solutions during 
the matter-dominated era, i.e at early times. Therefore, we consider p^, <C p m and p a , \p„\ <§; p m , or equivalently 
X0 <C 1 and \x a \ <C 1. Under these approximations, equations (|25| . ([26|) are simplified as: 



y 



W (0 ) 3 d^(^) 



(0) 



2 // 

a a (0) 



baa'. 



(0) 



3 dV^tr) 







= 0. 



(27) 
(28) 
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The zero subscript in parentheses "(0)" denotes just these zeroth-order solutions in terms of x<j,,Xcr, or equivalently 
in terms of the quintessence and phantom energy density, and must not be confused later on with the subscript 
without parentheses which stands for the present value of a quantity. 

Equation ([77]) can be easily solved analytically in the case of power-law potential (flU)) . The general solution for 
a<j> 7^ (which is the trivial case of a constant potential) is 
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(o) = C (a )a 2 +^, (29) 
where the function C^(a^,) is related to the potential parameters through 

r f n \ r 2a (2 + q ) 2 K ] 
L 3 / o„ l0 (4 + a-0) J 

In (|29| we have kept only the solution part that remains small (together with its derivative) for small a's, in order 
to be consistent with the matter-dominated approximation (p^Jp^l <C p m ). This is the reason for the absence of 
initial-condition dependent constants, which is just the central idea of tracker solutions. In addition, it is easy to see 
that this solution remains regular for — > —2. The solution for = — 2 can be shown to diverge at small a and 
thus we do not write it explicitly. 

Similarly, the general solution of equation (j2"8)) under potential (HU) is 

o"(0) = Ca(a a )a^, (31) 

where 

C a {a a ) = — — — . (32) 

L 3/9 m0 (4 + a CT ) J 

The zeroth-order solution for p^ can be calculated from ([!§]) under <C 1, that is p<p(o) = x <f>(o)Pm + ^('/'(o)): 
where x^q) — (acj)'^) 2 /6. The result is: 



[C<t>(a<f>)] 
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P0(o)=3p mO — 7T- rfl 2+a f- (33) 

«0(2 + a^j 

Furthermore, using the zeroth-order approximation of (|20|). i.e p CT (o) = ^ CT (o)Pm + ^((7(0)) with x CT ( ) = ~( acr (o)) /6> 
we obtain 

Mo) = -3ftnO Qff(2 + a|r) a — • (34) 

Similarly to the simple phantom case [l3| , expression (|34j) and the requirement of phantom-energy positivity leads 
to the constraint —2 < a CT < 0, that is the phantom potential is a normal power law. Additionally, (|33|) leads to 
< — 2 or > 0, with the second case being physically more interesting since it corresponds to the well-studied 
inverse power-law potential of the literature [H, [ill . [l2| . 

Finally, using the zeroth-order approximation of (f2"Tj) and (|22[) . namely p^(o) = £</>(o)Pm — ^(^o)) an d _p CT (o) = 
x a(o)Pm — Kt(< t (o)) respectively, we obtain 

; [QK)] 2 

'a^(2 + a^) 2 

P<t(o) = 6p m0 — 7^- 72 a 2+Q " ■ 36 



P^(o) = -6PmO ^ , _,_ ^ , 2 a 2+Q * (35) 



At this stage, the zeroth-order, that is early-time, results coincide with those of simple quintessence [12j or simple 
phantom models (l3j |. This was expected since the exact equations of motion (j2"5j) . (f2"6")l . which correspond to a coupled 
system for the two fields, under the zeroth-order approximation lead to the un-coupled equation system ([27|) . (|28p . 
Thus, at early times the two fields evolve independently. However, even at this approximation level, the complexity 
of the quintom model leads to qualitatively different results. Indeed, dark energy is attributed to both fields under 
definition (|12p . Therefore, we acquire: 



P<P{0)+Pa(0) Z 0^(2+^)2 " Z Q(t (2+q„)2 " 

; (0) = = 3^ ■ \"'J 

P0(o)+p CT(o) [c^W \c„M] 2 n -^- 



As we observe, although in the total absence of one of the fields this expression gives the corresponding results of 



simple quintessence (io™ 



|12| ) and simple phantom (i 



2+a+ M^i; — r~ t"«~- — v^(0) - 2 + a „ 

and despite the fact that they are not coupled at this zeroth-level, the behavior of W(o) is qualitative different. In 
particular, as can be clearly seen, is not a constant but it is a-dependent although it is the tracker solution. On 
the contrary, in simple quintessence and simple phantom scenarios one should go beyond tracker solutions, that is to 
intermediate times, in order to reveal a varying equation-of-state parameter. This feature makes the model at hand 
cosmologically interesting. We mention that according to the specific model and choice of parameters, u>(o) can evolve 
below or above —1, or even experience the phantom-divide crossing. 

Finally, we can calculate the zeroth-order behavior of the dark-energy density parameter as: 



[13]), when both fields are present 



n 



DE{0) 



P4>(o) + At(o) _ „ [C^Q^)] 



= 3 



3 



(38) 



where we have used that p m = p m o / a 3 . 

Expressions (|29|) -(|38 |) are the tracker solutions for quintom cosmology in power- law potentials. Equivalently they 
can be expressed as a function of time, considering a cx t 2 ! 3 since we are in the matter-dominated era, or as a function 



of the rcdshift z through ^ = 1 + z, with ao 
behavior of the quintom scenario as long as p, 



■ 1 the present value. They provide an excellent approximation to the 
\Pt/>\ Pm and Pa, \Pa\ "C Pm, th a t is at early times. 



III. COSMOLOGICAL EVOLUTION AT INTERMEDIATE TIMES 

As time passes and cosmological evolution continues, the scalar fields increase and dark-energy density becomes 
non-negligible, although still dominated by the dark-matter one. Therefore, we expect that progressively the various 
quantities will start diverging from the expressions obtained above. We are interested in the analytical investigation 
of this intermediate evolution stage, and thus we perform a first-order perturbation to the zeroth-order solutions of 
section |nl The obtained solutions will be significant for the period during the transition from matter to dark-energy 
domination. The corresponding quantities are denoting by the subscript "(1)", and the total ones by tilde, i.e: 



= 0(0) + 0(1), f = O-(O) + O-(l). P4> = P<t>(0) + P$(l), Pv = Pa(0) + P<j(l), W = W (0) + W {1 y (39) 

The aforementioned perturbations can be equivalently considered as keeping terms up to first order in x^, and x a (i.e 
keeping only awo) an d ^(o)) an d i n 0(i)/0(o) an d <T (i)/ cr (o)j i n the exact evolution equations (|25jl . ()26|) . An additional 
convenient formula can arise form the expansions of the potentials as V^(4>) — V^(4>^) + ~^~~' 0(i) 



similarly for V a (a). Substituting the expansions (f3"9")) into 
at this intermediate-time region: 



■0(0(1)), and 
we acquire the approximate evolution equations 



3k 



0(0) 



+ 



^(0)]^(0)}- 



(1)" 



■ 2 ( 5f7 (l) ~~ 3 I X 0(O) + X a(0)W(Q) } - 



3o^ ( o)[V*(0) + V a (a)} 3[p m + Pa(0) ] dV^ (0) ) 3 d 2 V^ (0) ) 

. ( 

2pm Pm 2 d<t>iQ) p m dej) 2 ^ 

3a^ ( o)[W) + V a {a)] 3{p m + p a(0) ] dV a (a {Q) ) 3 d 2 V a (a(o)] 



2p n 



Pm 



da 



where x 



0(0) 



finally obtain 



(0) 



daf Q) 



(i) = 0: 
(40) 

o-(i) = 0, 



(°0'(O)) 2 / 6 aild X cr 



(0) 



(41) 

(acr^) 2 /6. Inserting the explicit power-law potential forms ([13]) and (fT4| we 



and 



50(i) - v 



'(0)^(0) 



2(2 + a^pmo 
3a 4 



2pm0 



— a CT // 



(0) 



PrrM 12 

6a a (0) 



n 2 

PmO 



3a^(l + a^a 3 ^ 



-2-a„ 



■'(0) 



0(D 



PmO 



(42) 



'(i) 



5ct n + 



2 



a 

T 



'(O) ^(0) 



2(2 + a a )pma 
3a 4 



2pm0 



(0) a (o) 



£m0 -A/2 
6a ^(0) 



n 2 

PmO 



3av(l + a a )a 3 K a a 



(0) 



-X — Oto 

(0) 



PmO 



(43) 
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Despite their complicated form, these equation allow for a relatively simple general solution, if we keep only the part 
that remains small (together with its derivative) for small a's, in order to be consistent with the matter-dominated 
approximation. For > and —2 < a a < it reads: 



<j) (1) = B^a^a,,) a 2+ °*, (44) 

9 

<7(i) = B (T (a ( f > , oto) a 2 + a " , (45) 

where the functions B^^a^, a a ) and B cr {a l f )1 a a ) are related to the potential parameters through 

D/ \ o^-w \ {-[C^{a^)] 2 a a (2 + a CT ) 2 (6 + au) + [C a {a a )] 2 a^{2 + au)[12 + 4(a0 + a CT ) + a^}} 

B^OLf, Q. a ) = dC^a^) — — ■ — r ( 4b ) 

a^(2 + a^)a <T {2 + a<j) (28 + 8a^ + a^) 

R , x „~ , > {-[q(a CT )] 2 a (2 + a ) 2 (6 + a g ) + [C (a )] 2 « g (2 + a g )[12 + 4(a + a a ) + a cy]} 

^(a^.cv) = -6C a (a a ) — — — ■ — — , (47) 

a a {2 + a a )a^{2 + a^(28 + 8a a + a%) 

(that is B a (a,p,a a ) arises form B^(a,j,,a a ) changing the index <fi to a and vice versa, and adding an overall minus 
sign). In these expressions the complexity of the quintom model becomes clear, since they are qualitatively different 
from those of simple quintessence jl2| or simple phantom [l3T ] models. Indeed, since the first-order approximated 
evolution equations (|4"2"| . (|4"5|) are coupled, their general solution reveals a strong interdependence of the two fields, 
too. This feature will become even stronger in the w^-solution below, which by definition it already depends on both 
fields. Finally, in the total absence of one field, expressions (|46|) , (|47|) give exactly those of the corresponding simple 
models [IlOj]. 

The perturbation of p^, can be calculated from (|19p keeping the corresponding terms, thus: 

( f > (i) 

P<I>(1) = £0(O)P0(O) + £0(l)Pm ~ Q!0-7 V^,{4>), (48) 

0(0) 

where £0(i) = g a2< />(o)0(i) as ^ arises from x^ = x^o) + 2^(1) = a2 (0( 2 ) + ^0(o)0(i))/^ - Therefore, using also (|29[) . 
([55)) and (UH) in order to express the result in terms of Pmq), we obtain 

r 3[C (a0)] 3 -B0(a ,a (T )Q !0 (a 2 -4)-, b_ 

= 1 2C,K)(2 + ^) 2 H 0) ° * " (49) 

Repeating the same procedure for p a we acquire: 

f 3[C a (a a )} 3 + ^(a , a tT )a CT (a 2 - 4) ^ e 
Pff(1) = "1 2a(a g )(2 + QCT ) 2 K<°> a2+ ^ • (50) 

The perturbation of w can arise form (|12[) as 

_ P0(1) +Pa(l) P<f,(0) + P«r(0) r , , , 

; T 1 12 ^(1) + Pct(i) j I 51 ) 

/V(0)+P<x(0) W(0) + Pcr(0)\ 

with p^(i) = X0( O )P0( O ) + x<p(i)p m + Q; 0^■^( ( / , ) ( and similarly for p CT (i)). The final result reads 

a 2 (2 + a 2 )a 2 (2 + a 2 )a 

W(i) = — ^ -^(di+d 2 + d 3 + d 4 ), (52) 



where 



2 [C CT (a CT )] 2 a (2 + c^)a 2+c v - [C (a )] 2 a CT (2 + a CT ) a — 



[C (q )] 3 (4 + a ){3[C0(a )] 3 + ^(a^, a CT )a0(2 + a )(6 + a )} 

a 2 ^(2 + a ) 5 



, [C (T (a (T )] 3 (4 + a (T ){-3[C' (T (a (T )] 3 + J B (T (a ,a (T )a (T (2 + a (T )(6 + a (T )} eg-^) 

do — — — a 2 + q ct 

a 2 (2 + a ff ) 5 



rf 3 
e?4 



CV(o^)[C g (a g )] 2 {3(4 + a tT )[C (a )] 3 + ff (a , a CT )a (2 + a )[24 + (10 + a<f> )a a }} ^l. 

«0(2 + a (j) ) 3 a (T (2 + a a ) 2 
C CT (a (T )[C^(a )] 2 {™3(4 + a )[C CT (a (T )] 3 + S CT (a , a CT )a CT (2 + a a )[24 + (10 + a CT )a ]} 



,2+a A 2+ a 



a a {2 + a CT ) 3 a0(2 + a^) 2 
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Note that expressions (|49|) , (f50|> and (|52|) in the total absence of one field, coincide with those of the corresponding 
simple quintessence |12| and simple phantom [l3| models. The qualitative difference of quintom scenario is also 
reflected in the un-specified sign of the correction terms, contrary to simple models where these corrections have a 
constant sign (</>m is always negative in both quintessence and phantom cases, and p<f>(i) > and ^^(i) < in simple 
quintessence while p^i) < and u>0(i) > in simple phantom scenario). 

In summary, inserting ([44]) . (|45|) . (|49|) . (|50|) and ([52]) in the perturbative expansion (l39|) we finally obtain the 
analytical solutions for intermediate times, namely <f>, cf, p^, po- and u;. Finally, fius can be calculated as £Ide = 
(P4> + Pcr)/(p0 + /V + Pm)- We mention that since we have obtained the aforementioned quantities as a function of 
the scale factor, it is straightforward to express them as a function of the redshift z through - = 1 + z. 



IV. COMPARING ANALYTICAL AND NUMERICAL RESULTS 



In the previous two sections we extracted analytical expressions for the quintom scenario in general power-law 
potentials, when the universe is dominated by the dark matter sector. Concerning their applicability, we argued that 
these solutions are valid at early and intermediate times, that is at both high and low redshift. In order to examine 
the precision of our results, we evolve numerically the exact cosmological system calculating the exact quantities </>(a), 
Qde(o-), w(a), and then we study their divergence from the zeroth and first order formulae derived above. This can be 

achieved by examining the ratios <fi(o)/4>i 4>I4>, cr (o)/ <J : and similarly £Ide(o)/^de, ^de/^de and w^/w, w/w. 
The initial conditions are chosen consistently with initial matter domination, since if this requirement is fulfilled then 
the results do not depend on the specific initial conditions, as expected. Finally, we fix and K a in order to acquire 
fl m o ~ 0.28 and £Ideo ~ 0.72 at present. Note that contrary to simple quintessence and simple phantom cases, one 
can satisfy these observational requirements by many different realizations of the quintom scenario (i.e many different 
potential-parameter choices), which is an additional advantage of the model at hand. 

In fig. [T]we present the ratios of the zeroth and first order field solutions to the exact numerical values, as a function 
of the redshift. The calculations have been performed for two potential combinations, namely = 2, a a = —I and 




FIG. 1: (Color Online) Upper graph: The ratios 4>^/<j) (upper two curves) and <j>/4> (lower two curves) as a junction of the 
redshift, for — 2,a a = —1 (black, solid) and aj, = l,cta = —0.5 (red, dashed). Lower graph: The ratios er^/a (upper two 
curves) and a /a (lower two curves) as a function of the redshift, for = 2,a a — —1 (black, solid) and ct^ = l,a a = —0.5 
(red, dashed). 



= l,a a = —0.5. As we observe, <j>iQ)/<p and cr( )/cr are very close to 1 for z > 1.5, that is the tracker solution is a 
very good approximation at this early evolution stage, with errors of less than 3%. Furthermore, for 1.5 > z > 0.5 
the zeroth-order solution is not a good approximation (with error 15%) but the first-order one remains within 95% 
accuracy. As expected these errors are larger than the corresponding results of simple quintessence [l2l ] and simple 
phantom [13] models, since the exact system is coupled and thus the evolution depends on both fields. Finally, at 
late times, that is for z < 0.5, the scalar fields are enhanced significantly, dark energy dominates the universe, and 
our approximation breaks down rapidly. 

In fig. [5] we depict ^de(o)/^de and Q,oe/^de as a function of z. The divergence of the zeroth-order solution 
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FIG. 2: (Color Online)27ie ratios Qde(o)/Qde (upper two curves) and Qde /Qde (lower two curves) as a function of the 
redshift, for otj, = 2, a CT = —1 (black, solid) and = l,a CT = —0.5 ('red, dashed). 



from the exact one starts at z w 1.5. However, we can see that the first-order solution is very accurate (with error 
less than 2%) up to z m 0.5. After that point, the first-order solution starts diverging rapidly from the exact one and 
our approximation is not valid. 

Finally, in fig. [3]we present w^/w and w/w as a function of the redshift. We observe that the tracker solution is 
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FIG. 3: (Color Online) The ratios W(p)/w (upper two curves) and w/w (lower two curves) as a function of the redshift, for 
— 2, ct a = —1 (black, solid) and = 1, a ff = —0.5 (red, dashed). 



very accurate up to z ~ 1.5, while the first-order one agrees with the exact evolution within 2% up to z ~ 0.5. As 
expected these errors are larger than the corresponding results of simple models [TH, [l3| , since the error of w depends 
on the errors of both fields. 

Therefore, from these tree figures we conclude that the tracker, that is the zeroth-order, solution is accurate within 
an error of 2% at early times and up to z w 1.5. At intermediate times, that is at low redshifts, we have to use the 
first-order analytical solution, which agrees with the exact cosmological evolution up to z as 0.5, within an error of 
2%. After that point, the dark energy sector enhances, it dominates the dark matter one, and our approximation 
breaks down rapidly. 
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V. COSMOLOGICAL IMPLICATIONS AND DISCUSSION 

Having tested the accuracy of our analytical expressions and having determined their applicability region, we can 
use them to describe an arbitrary quintom evolution in power-law potentials up to z — 0.5. Clearly, the tracker 
solutions are much simpler and can be used as a first approximation. Contrary to simple quintessence and simple 
phantom models, w presents a varying behavior even at this zeroth-order approximation level, and thus it can offer 
a reliable description of the system (quantitatively up to 2 « 1.5 and qualitatively later on). For intermediate times 
one can use at will the analytical expressions for w and 1]^, which are very accurate up to z ~ 0.5. Here we present 
an additional simplification which allows for an easier application. Observing that their complicated form is due to 
the a-independent coefficients (which are determined by the potential parameters) rather than the exponents of a, 
we can make an easy shortening keeping the most significant a-terms and thus resulting to a convenient simulation 
of w-evolution. Therefore, expressing the formula in terms of the redshift z through i = 1 + z, we finally acquire: 

w fltz (z) = u>i(l + z) Sl + (w - Wl )(l + z) s \ (53) 

In this w-parametrization wq is its present value, which can be taken from observations, while w\, s\ and S2 are 
determined by the potential exponents in a rather complicated way. Equivalently, we can consider w±, si and S2 
as additional, free parameters. In this way, the relative large number of parameters in the fit (|53[) is attributed to 
the general behavior of quintom scenario, which must be able to present a large class of cosmological evolutions. 
Furthermore, the use of three or four parameters in w-parametrizations is usual in the literature [Til ]. 

Let us examine the accuracy of the fit (|53[) . comparing it with the exact cosmological evolution of the model at 
hand. In fig. 3]we depict the ratio Wfuz/w as a function of z for two potential choices. In this figure the value wq has 
been determined numerically, that is why Wfit z and w coincide both initially and finally. As we observe, the proposed 
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Z 

FIG. 4: (Color Online) The ratio Wfuz/iv, where Wfitz is given by \53\) . as a function of z, for = 2,a a — —1 (black, solid) 
and a?0 = l,a CT = —0.5 (red, dashed). 

w-parametrization, given by (|53p . agrees with the exact evolution within 99%. Similarly to the figures of the previous 
section, the errors are larger for larger a§ and \a a \, however small a^, \a a \ are expected to be more realistic. The 
very small error is due to the general and multi-parametric form of the fit (|53p . which allows for a very satisfactory 
fit of a general quintom evolution in power-law potentials. 

Another useful approximated relation, that can arise from our analytical results, is the expression of w as a function 
of CIde- Indeed, since we have derived the relations w(a) and £Ide{o-), we can keep the most significant a-terms and 
then eliminate a, resulting to w(Cl de)- Doing so we obtain: 

WfttQ.DE (^De) = W 2 + W 3 ^DE + W ^ 8 DE- ( 54 ) 

The various parameters can be determined by the specific potential choice in a rather complicated way. Equivalently 
we can consider them as free parameters, suitably fitted in order to acquire the best agreement with the exact 
cosmological evolution. Using a \ 2 minimization routine to determine the best-fit values of these parameters, in fig. 
[S]we present the ratio Wfitn DE l w as a function of flnE for two potential choices, where the present r^D^-value has 
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FIG. 5: (Color Online) The ratio Wfitci DE /w, where Wfitn DE is given by (5^[ j, as a function ofSloE, for a.^ = 2, ( 
(black, solid) and = 1, <x a = —0.5 ('red, dashed). 



been taken to be 0.72, both in Wfitn DE an d w - The error introduced by the fit (|54[) is less than 2.5% for the examined 
potential cases, mainly at late times. The fact that it is larger than the corresponding error of Wfi tz , is a result of 
the simplifications we performed in order to eliminate a and obtain Wf anDE^DE)- Finally, we stress that fits (|53p 
and (fS"i|) can be considered to be valid up to present epoch (z = and FIde ~ 0.72) although they have arisen from 
expressions valid up to z ss 0.5, since the relative large number of free parameters leads to small errors and one can 
always adjust the fit to coincide with the exact evolution both initially and finally. 

Since we have examined the accuracy of our w(z) and w(Qde) parametrizations, we can use them at will to 
describe a large class of cosmological evolutions. In fig. [5] we depict Wfu z (z) for four different cosmological cases. 
Case I corresponds to a quintessence-like evolution, where w remains always larger than —1. Case II corresponds to 
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FIG. 6: (Color Online) Wf itz (z) for four different cosmological cases, as it arises from parametrization S53\) . Case I: wo = —0.9, 
wi — —1, si = —1/3, S2 = 1/5 (black, solid), case II: wo = —1.05, wi — 1, si — —1/3, S2 = 1/5 (red, dashed), case III: 
wo = —1.05, wi — —1, si = —1/3, S2 = 1/5 (green, dotted), case IV: wo = —1.05, wi — —1.8, si = 1, S2 = 1.8 (blue, 
dashed- dotted). 



a phantom- like evolution, with w < — 1. Case III corresponds to a quintom evolution where w crosses the phantom 
divide in the recent cosmological past, as might be suggested by observations within 95% confidence level [J. It is 
interesting to see that (|5T{|) allows for a non-monotonic Wfi tz (z), too. Indeed, case IV corresponds to such a behavior, 
where w experiences the — 1-crossing from above to below and then it starts increasing again. 

In summary, as was mentioned above, expression (|53p can describe a very large class of cosmological evolutions. 
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Similar results can be obtained using w/itn DE (^de) from (|54p . where one can acquire various w(flnE) cosmological 
behaviors. Finally, it is easy to see, that the case of simple quintessence can be described by ([55|) and (j5~4"|) using 

Sl = 0, S 2 = -24^, «3 = 1, tfll = t02 = -24^, ™3 - - (2+4x^+8^+28) ' WA = ' wMe f ° r the Sim P k P hantom 

model the corresponding parameters are the same with — > cvo-. However, it is more convenient to consider them as 
free parameters. 

VI. CONCLUSIONS 

In this work we studied quintom models with power-law potentials. Firstly, we extracted the tracker solutions 
under the assumption of matter domination. These solutions correspond to the general and common behavior of all 
such models at early times, that is at high redshifts. Contrary to simple quintessence and simple phantom cases, even 
at this zeroth-order approximation level, w is not constant but evolving, which is a result of the complex nature of 
quintom scenario. The tracker solutions are quantitatively valid up to z k 1.5 (within 98% accuracy), but since this 
zeroth-order approximation possesses the qualitative features of the exact solution (varying and not constant w), it 
could be used up to low redshifts in order to provide a first picture of the corresponding cosmological evolution. 

In addition, we extracted the general cosmological solutions at intermediate times, that is at low redshifts, which 
is the period during the transition from matter to dark-energy domination. Such a solution can be very useful in 
dark energy observations, since probes based on supernovae la, WMAP and SDSS, are related to this cosmological 
epoch [l[. The comparison with the exact evolution shows that these first-order solutions are accurate within 2% up 
toz« 0.5. 

In order to use the aforementioned results in a relative simple way, we have further approximated them extracting 
two new w-parametrizations, one as a function of the redshift z and one as a function of Qde- Although the 
fitting parameters can be expressed as a function of the potential ones, one can equivalently consider them as free 
parameters suitably chosen to describe an arbitrary evolution. Due to the relative large number of free parameters, 
such parametrizations are very accurate up to present epoch. Thus one can use them in order to describe various 
quintom evolution sub-classes, including quintessence- like or phantom-like cases, realization of the — 1-crossing, non- 
monotonic w{z) evolution etc. 

The above analysis shows that a two-field (one canonical and one phantom) quintom scenario with power-law 
potentials can offer a good description of dark energy evolution, in agreement with observations [lj. Moreover, the 
fact that power-law potentials can be justified through supersymmetric considerations [10] is an additional advantage. 
However, as it is usual in models where phantom fields are present, the quantum behavior of the examined scenario 
needs some caution, since the discussion about the construction of quantum field theory of phantoms is still open in 
the literature. For instance in fl5j the authors reveal the causality and stability problems and the possible spontaneous 
breakdown of the vacuum into phantoms and conventional particles, but on the other hand in [16| the phantom fields 
arise as an effective description, consistently with the basic requirements of quantum field theory. The comparison 
with observations and a robust theoretical justification should be the decisive tests for the present model. 



[1] A.G. Riess et al. [Supernova Search Team Collaboration], Astron. J. 116, 1009 (1998); S. Perlmutter et al. [Supernova 
Cosmology Project Collaboration], Astrophys. J. 517, 565 (1999); A. Melchiorri et al, Astrophys. J. Lett. 536, L63 (2000); 
D. N. Spergel et al, Astrophys. J. Suppl. 148, 175 (2003); Bridle, O. Lahab, J. P. Ostriker and P. J. Steinhardt, Science 
299, 1532 (2003); J. L. Tonry et al, Astrophys. J. 594, 1 (2003); S. W. Allen, et al., Mon. Not. Roy. Astron. Soc. 353, 
457 (2004). 

[2] P. Binetruy, C. Deffayet, D. Langlois, Nucl. Phys. B 565, 269 (2000); G.R. Dvali, G. Gabadadze, M. Porrati, Phys. Lett. B 
485, 208 (2000); S. Capozziello, Int. J. Mod. Phys. D 11, 483 (2002); S.Nojiri and S. D. Odintsov, Phys. Rev. D 68, 123512 
(2003); A. Lue and G. Starkman, Phys. Rev. D 67, 064002 (2003); P. S. Apostolopoulos, N. Brouzakis, E. N. Saridakis 
and N. Tetradis, Phys. Rev. D 72, 044013 (2005); G. Calcagni, S. Tsujikawa and M. Sami, Class. Quant. Grav. 22, 3977 
(2005); S.Nojiri and S. D. Odintsov, Int. J. Geom. Meth. Mod. Phys. 4, 115 (2007); R. P. Woodard, Lect. Notes Phys. 
720, 403 (2007); F. K. Diakonos and E. N. Saridakis, JCAP 0902, 030 (2009). 

[3] P. J. E. Peebles and B. Ratra, Astrophys. J. 325, L17 (1988); B. Ratra and P. J. E. Peebles, Phys. Rev. D 37, 3406 (1988). 

[4] C. Wetterich, Nucl. Phys. B 302, 668 (1988); M. S. Turner and M. White, Phys. Rev. D 56, 4439 (1997); R. R. Caldwell, 
R. Dave and P. J. Steinhardt, Phys. Rev. Lett. 80, 1582 (1998); Z. K. Guo, N. Ohta and Y. Z. Zhang, Mod. Phys. Lett. A 
22, 883 (2007); O. Hrycyna and M. Szydlowski, Phys. Rev. D 76, 123510 (2007); S. Dutta, E. N. Saridakis, R. J. Scherrer, 
larXiv:0903.34i"2l [astro-ph.CO]. 

[5] I. Zlatev, L. Wang, and P. J. Steinhardt, Phys. Rev. Lett. 82, 896 (1999); A. R. Liddle and R. J. Scherrer, Phys. Rev. D 
59, 023509 (1999); P. J. Steinhardt, L. Wang, and I. Zlatev, Phys. Rev. D 59, 123504 (1999). 



12 



[6] R. R. Caldwell, Phys. Lett. B 545, 23 (2002); R. R. Caldwell, M. Kamionkowski and N. N. Weinberg, Phys. Rev. Lett. 91, 
071301 (2003); V. K. Onemli and R. P. Woodard, Phys. Rev. D 70, 107301 (2004); P. F. Gonzalez-Diaz and C. L. Siguenza, 
Nucl. Phys. B 697, 363 (2004); M. Sami, A. Toporensky, P. V. Tretjakov and S. Tsujikawa, Phys. Lett. B 619, 193 (2005); 
A. Vikman, Phys. Rev. D 71, 023515 (2005); Z. K. Guo, R. G. Cai and Y. Z. Zhang, JCAP 0505, 002 (2005); S. Nojiri and 
S. D. Odintsov, Phys. Rev. D 72, 023003 (2005); H. Garcia-Compean, G. Garcia- Jimenez, O. Obregon, and C. Ramirez, 
JCAP 0807, 016 (2008); E. N. Saridakis, larXiv:0811.1333l [hep-th ]; M. Szydlowski and O. Hrycyna, JCAP 0901, 039 
(2009); X. m. Chen, Y. Gong and E. N. Saridakis. larXiv:0812. 11171 [gr-qc], 
[7] B. Feng, X. L. Wang and X. M. Zhang, Phys. Lett. B 607, 35 (2005); Z. K. Guo, et ai, Phys. Lett. B 608, 177 (2005); 
H. Wei, R. G. Cai and D. F. Zeng, Class. Quant. Grav. 22, 3189 (2005); M.-Z Li, B. Feng, X.-M Zhang, JCAP, 0512, 
002 (2005); P. x. Wu and H. w. Yu, Int. J. Mod. Phys. D 14, 1873 (2005); R. Lazkoz and G. Leon, Phys. Lett. B 638, 
303 (2006); Y. f. Cai, H. Li, Y. S. Piao and X. m. Zhang, Phys. Lett. B 646, 141 (2007); W. Zhao and Y. Zhang, Phys. 
Rev. D 73, 123509 (2006); B. Feng, M. Li, Y.-S. Piao and X. Zhang, Phys. Lett. B 634, 101 (2006); M. R. Setare and 
E. N. Saridakis, Phys. Lett. B 668, 177 (2008); M. R. Setare and E. N. Saridakis, JCAP 0809, 026 (2008); M. R. Setare 
and E. N. Saridakis, Phys. Lett. B 671, 331 (2009). 
[8] P. S. Apostolopoulos, and N. Tetradis, Phys. Rev. D 74, 064021 (2006); H.-S. Zhang, and Z.-H. Zhu, Phys. Rev. D 
75, 023510 (2007) ; Y. Gong an d A. Wang, Phys. Lett. B 652, 63 (2007); E. N. Saridakis, P. F. Gonzalez-Diaz and 
C. L. Siguenza. [arXiv:0901. 12131 [astro- ph], 
[9] R. J. Scherrer and A. A. Sen, Phys. Rev. D 77, 083515 (2008); R. J. Scherrer and A. A. Sen, Phys. Rev. D 78, 067303 
(2008); M. R. Setare and E. N. Saridakis, Phys. Rev. D 79, 043005 (2009). 

[10] P. Binetruy, Phys. Rev. D 60, 063502 (1999); A. Masiero, M. Pietroni and F. Rosati, Phys. Rev. D 61, 023504 (2000). 

[11] J. P. Kneller and L. E. Strigari, Phys. Rev. D 68, 083517 (2003); L. R. W. Abramo and F. Finelli, Phys. Lett. B 575, 
165 (2003); X. Zhang, Mod. Phys. Lett. A 20, 2575 (2005); M. Yashar, B. Bozek, A. Abrahamse, A. J. Albrecht and 
M. Barnard. larXiv:0811. 22531 [astro-ph] . 

[12] C. R. Watson a nd R. J. Scherrer, Phys. Rev. D 68, 123524 (2003). 

[13] E. N. Saridakis. la"rXiv:0902. 39781 [gr-qc] . 

[14] G. Efstathiou, Mon. Not. R. Astron. Soc. 310, 842 (1999); A. R. Cooray and D. Huterer, Astrophys. J. 513, L95 (1999); 

M. Chevallier and D. Polarski, Int. J. Mod. Phys. D 10, 213 (2001); E. V. Linder, Phys. Rev. Lett. 90, 091301 (2003). 

E. Di Pietro and J. F. Claeskens, Mon. Not. Roy. Astron. Soc. 341, 1299 (2003); M. Szydlowski, A. Kurek and A. Krawiec, 

Phys. Lett. B 642, 171 (2006); E. J. Copeland, M. Sami and S. Tsujikawa, Int. J. Mod. Phys. D 15, 1753 (2006); H. Wei 

and S. N. Zhang, Phys. Lett. B 644, 7 (2007). 
[15] J. M. Cline, S. Jeon and G. D. Moore, Phys. Rev. D 70, 043543 (2004). 

[16] S. Nojiri and S. D. Odintsov, Phys. Lett. B 562, 147 (2003); S. Nojiri and S. D. Odintsov, Phys. Lett. B 571, 1 (2003). 



